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ABSTRACT 

A new aspect of the vacuum structure of 2+1- dimensional Thirring model is presented. 
Using the Fierz identity, we split the current- current four-Fermi interaction in terms of a 
matrix valued auxiliary scalar field and compute its effective potential. Energy consideration 
shows that contrary to earlier expectations, parity in general is spontaneously broken at any 
finite order of N, where N is the number of the two component spinors. In the large N limit, 
there does not exist a stable vacuum of the theory thereby making the application of the 
large N limit to Thirring model dangerous. A detailed analysis for parity breaking solutions 
in N=2,3 cases is given. 
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It is well known that 2+1-dimensional QED admits dynamical mass generations. [1]- [3] 
Solving the Dyson-Schwinger gap equation of the theory, Appelquist et al. [4] have argued 
that in the large N limit, where N is the number of two-component fermions, masses are 
dynamically generated in such a way to preserve the overall parity even though each indi- 
vidual mass term violates parity symmetry. Recently, similar analysis has been applied to 
2+1-dimensional Thirring model of vector- vector type four- Fermi interactions with a same 
type of dynamical mass generations. [5] The four- Fermi interaction model in 2+1 dimensions 
are also known to be renormalizable in the 1/N-expansion. [6] 

In this Letter, we bring attention to the danger of this type of analysis, particularly 
in the case of 2+1-dimensional Thirring model. We analyze symmetry-breaking patterns 
in 2+1-dimensional Thirring model for finite N as well as for the large N limit. Using 
the Fierz identity, we split the current- current four-Fermi interaction in terms of a matrix 
valued auxiliary scalar field and instead of solving the Dyson-Schwinger gap equation, we 
compute the effective potential of the auxiliary field to study the dynamical mass generations. 
Energy consideration shows that contrary to earlier expectations, overall parity in general is 
spontaneously broken at any finite order of N. In the large N limit, there does not even exist 
a stable vacuum of the theory which brings doubts to the 1/N-analysis of 2+1-dimensional 
Thirring model. The discrepancy of our result with earlier works [4] based on the Dyson- 
Schwinger equation arises because in [4], the overall parity-breaking solution was discarded 
in the large N limit analysis since it exceeded the ultraviolet cut-off. However, our result 
shows that for finite N and for certain range of coupling constants the parity-breaking solution 
becomes a true vacuum. In the large N limit, it exceeds the ultraviolet cut-off as before which 
however makes other solutions semi-classically unstable. We argue that similar difficulty also 
arises for 2+1-dimensional QED. 

Consider the 2+1-dimensional Thirring model given by the Lagrangian 



where ip 1 are two-component spinors and i runs over from 1 to N. In two-component repre- 
sentation, the Dirac 7 matrices are given in terms of the Pauli matrices, 



(1) 



7 = ct 2 , 7 1 = i&3 , 7 2 = «cti • 



(2) 
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Note that the Lagrangian in Eq.(l) is invariant under the global U(N) transformation as 
well as under the parity transformation P, 



P:(x,y,t) -> (-x,y,t) 

if>'(x) -»• i/>(x) = P^{x)P' 1 = a^{x) . (3) 

It is believed that due to the strong coupling of fermions, the parity symmetry breaks 
down spontaneously at the quantum level which generates dynamical masses for fermions. 
However, it was argued that for even number of fermions iV = 2n, half of the fermions 
ip 1 ;i = 1, n get mass m while the other half ; % — n + 1, 2n get mass — m so as to 
preserve the overall parity P4 = PZ2 where Z 2 mixes fermions: [5] 

Z 2 ^(x) -> ^ n+i (x) , Z 2 ^ n+ \x) -> P(x) ; i = 1, n . (4) 

This argument is based on the Dyson-Schwinger equation method in the large N limit which 
gives rise to the same type of dynamical parity breaking for 2+1-dimensional QED.[4] In 
the following, we show that contrary to the above argument the overall parity P4 is in fact 
broken. To do so, we first note that the action Eq.(l) can be written via Fierz transformation 

L = pi-fdrf + ^( W + j^VW (5) 

where the four-Fermi interactions can be splitted by introducing a matrix-valued auxiliary 
field Mij 

L = + M)ip - ?-trM* + 4g{ * + 2) {trMf . (6) 

This reduces to the original Lagrangian when is eliminated by integration which identifies 
M,^ with Mij = j^i[) k il) k 5 1 ^ + ip 1 . In order to understand the struture of vacuum, we 
compute the effective potential for M^ which amounts to summing the one loop fermion 
diagrams at zero momentum. Since the matrix M is hermitian, M may be diagonlized with 
real components M = diag(\i, X 2 , Ajv) and the effective potential, invariant under the 
diagonalization, is given by 

(7) 
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where a cut-off has been introduced, after Wick rotation, at k 2 = A 2 in order that the 
integral be well defined. In terms of the dimensionless quantity x^ = \/A and the rescaled 
coupling constant g = 2Ac//37r 2 , Eq.(7) becomes 

N N AT N AT N 

Veff = - E Ml + *?) ~ 2 E \^\ 3 (tan^\x t \ - ) + ( - 2) £ *? - ^tt^t(E ^) 2 • (8) 



i=i i=i 



In the following, we analyze various symmetry-breaking patterns of the potential V e ff ac- 
cording to the value of coupling constant g. First consider the case N = 2. The classical 
stability of the perturbative vacuum (xi = 0, x 2 = 0) is governed by the second derivatives 
of the potential V e ff at (0, 0): 

a- 92 V, r 3 D _ d 2 V 1 ^_d 2 V. r 3 /n . 

A = m l(o - o) = " 6 + r 9^ l(o ' o) = "i ' c = M Vm + s ( } 

which imply that the perturbative vacuum is a local maximum if 5 2 — AC < and A+C < 0, 
a local minimum if B 2 — AC < and A + C > 0, or a saddle point if -B 2 — AC > 0. Thus, 
the potential V e ff(xi, X2) and the perturbative vacuum (0,0) have the following properties; 

i) g < 0: V e ff is bell-shaped and (0, 0) is the absolute maximum of V e ff so that there is 
no stable vacuum. 

ii) < g < 1/3: V e ff is cup-shaped and (0,0) is the absolute minimum of V e ff so that 
(0, 0) is a stable vacuum. 

hi) 1/3 <g< 2/3: There are two local minima at (m,m) and (— m, —m) for < m < 1. 
V e ff is double- well shaped when restricted on the line x\ = X2 and U-shaped when restricted 
on the line x± = —x<i and (0,0) is a saddle point of V e ff. Thus, (±m, ±m) become stable 
vacua which break parity symmetry spontaneously. 

iv) 2/3 < g: There are four local minima at (mi, mi), (—mi, —mi) and [m-i, —m-i), {—m-i, m^). 
V e ff is double- well shaped when restricted on the line x\ = X2 as well as on the line x\ = —x 2 . 
m i > m 2 and V e ff(±mi,±mi) < V e ff(±m 2 ,T m 2)- Thus (±mi,±mi) become stable vacua 
which break parity symmetry spontaneously. However numerical computation shows that 
for g > 1.5, mi > 1. This implies that the dynamically generated mass exceeds the cut-off 
A and therefore (±m 1 ,±m 1 ) are not sensible vacua. Nevertheless, this makes (±m 2 , ^2) 
semi-classically unstable so that there is no stable vacuum in this case. 

For N = 3, similar analysis leads to the following properites of the potential; 
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i) g < 0: (0, 0, 0) is the absolute maximum of V e ff so that there is no stable vacuum. 

ii) < g < 2/5: (0,0,0) is the absolute minimum of V e ff so that (0,0,0) is a stable 
vacuum. 

hi) 2/5 < g < 1: (0, 0, 0) is a saddle point and there are two local minima at (±m, ±m, ±m) 
which become stable vacua breaking parity spontaneously. 

iv) 1 < g: (0, 0, 0) is a local maximum and there are six local minima with same poten- 
tial value and two absolute minima at (±m, ±m, ±m) breaking parity spontaneously which 
become sensible vacua for g < 1.9. 

In order to understand the large N behavior, we first consider the even N case and then 
take the large N limit. For N = 2n, we look at the special sector of the domain of the 
potential: x± = x 2 = ■ ■ ■ = x n = x and x n+ i = x n+2 = ■ ■ ■ = x 2n = V where the symmetry 
breaking vacuua are expected to arise. Then the effective potential becomes 

^V eff (x, y) = -Ml + x 2 ) + ln(l + y 2 )] - 2[\x\ 3 (ta n - 1 \x\ - |) + ^{tan- 1 ^ - |)] 

+(f-^ 2 + . 2 )- (* + y?- do) 

In terms of a rescaled coupling constant g = 2g/N, the structure of the effective potential 
jjV e ff(x, y) is the same as that of the N = 2 case with the replacement of range of g : i) 
g < , ii) < g < 4/3(2 + N) , hi) 4/3(2 + N) < g < 2/3 , iv) 2/3 < g . For 
large N, the potential ^K// possesses local minima along the line x = — y when g > 2/3. 
However, in the limit N oo the potential becomes unbounded below when restricted 
on the line x = y for all values of g . Therefore, there does not exist a stable vacuum in 
the large N limit! At first sight, this result seems to be in contradiction to earlier works 
based on the Dyson-Schwinger equation method. In the Dyson-Schwinger type analysis 
with the 1/N- approximation, local minima arise along the line x = —y which agrees with 
our result. However another set of solutions to the Dyson-Schwinger equation, which arise 
on the line x = y and also in agreement with our results, were discarded because they 
exceed the momentum cut-off range. Our results show that these solutions can not be 
simply discarded since they make other local minima semi-classically unstable. Finally, we 
note that even though our analysis is only for the case of 2+1 dimensional Thirring model, 
similar vacuum instability in the large N limit may exist in other 2+1-dimensional model as 
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well. In particular the vacuum struture of 2+1-dimensional QED was essentially the same 
as that of the Thirring model in the Dyson-Schwinger approach which suggests that the 
1/N- approximation of 2+1-dimensional QED is as well dangerous. 
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